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Nonlinearity arising from noncooperative transcription factor binding 
enhances negative feedback and promotes genetic oscillations 

Ivan M. Lengyelp- Daniele SoroldoniAndrew C. Oates,® Luis G. Morelli * 1 2 3 * 


We study the effects of multiple binding sites in the promoter of a genetic oscillator. We 
evaluate the regulatory function of a promoter with multiple binding sites in the absence 
of cooperative binding, and consider different hypotheses for how the number of bound 
repressors affects transcription rate. Effective Hill exponents of the resulting regulatory 
functions reveal an increase in the nonlinearity of the feedback with the number of binding 
sites. We identify optimal configurations that maximize the nonlinearity of the feedback. 
We use a generic model of a biochemical oscillator to show that this increased nonlinearity 
is reflected in enhanced oscillations, with larger amplitudes over wider oscillatory ranges. 
Although the study is motivated by genetic oscillations in the zebrafish segmentation clock, 
our findings may reveal a general principle for gene regulation. 


I. Introduction 

Cells can generate temporal patterns of activity 
by means of genetic oscillations w- Genetic os¬ 
cillations are biochemical oscillations in the levels 
of gene products EHH]. They can be produced 
by negative feedback regulation of gene expression, 
in which a gene product inhibits its own produc¬ 
tion directly or indirectly [15]. Such autoinhibi¬ 
tion is often performed by transcriptional repres¬ 
sors, proteins or protein complexes that bind the 
promoter of a gene and inhibit the transcription of 
new mRNA molecules mu , see Fig. [TJ A theoret¬ 
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ical description of biochemical oscillations requires 
a delayed negative feedback, together with suffi¬ 
cient nonlinearity and a balance of the timescales 
of the different processes involved [T5] . Delays may 
occur naturally in transcriptional regulation, since 
the assembly of gene products involves intermedi¬ 
ate steps |16II19II20| . Nonlinearity refers to the pres¬ 
ence of nonlinear terms in the equations describing 
the dynamics. Such nonlinear terms may occur in 
the equations due to the presence of cooperative 
biochemical processes, where cooperativity is un¬ 
derstood as a phenomenon in which several com¬ 
ponents act together to orchestrate some collective 
behavior [21]. Although some processes giving rise 
to nonlinear terms are known, in general it is still 
an open question how nonlinearity is built into ge¬ 
netic oscillators. 

A compelling model system for genetic oscilla¬ 
tions is the vertebrate segmentation clock WI\W2 . 
This is a tissue-level pattern generator that controls 
the formation of vertebrate segments during em¬ 
bryonic development [25J[2E/ The spatiotemporal 
patterns generated by the segmentation clock are 
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thought to be initiated by a genetic oscillator at the 
single cell level 1711^1110112711^] . In this genetic oscil¬ 
lator, negative feedback is provided by genes of the 
Her family, which code proteins that form dimers 
and can act as transcriptional repressors [3(TH34j . 
The time taken from transcription, translation and 
splicing may introduce the necessary feedback de¬ 
lays in zebrafish and mouse Hama- One source 
of nonlinearity in the segmentation clock oscilla¬ 
tor is the dimerization of gene products that bind 
the promoter of cyclic genes [32]. However, this 
may be insufficient to generate the observed oscil¬ 
lations in the levels of gene products. One way to 
increase nonlinearity would be cooperative binding 
of repressors to regulatory binding sites at the pro¬ 
moter [2T[[38][39]- Cooperative binding to multi¬ 
ple binding sites can make the negative feedback 
steeper [50] • A similar effect occurs with ultra¬ 
sensitivity in phosphorylation cascades [HI . The 
presence of clusters of binding sites for transcrip¬ 
tion factors may be a common motif in gene regula¬ 
tion [42j . and cooperativity in transcription factor 
binding has been reported for some systems [551 . 

In zebrafish, multiple binding sites for Her dimers 
have been identified in the promoter region of 
Herl, Her7, and other genes of the Notch path¬ 
way [SHHSi]. However, there is no evidence for 



Figure 1: Delayed autoinhibition can produce ge¬ 
netic oscillations. (A) The gene (light blue box) 
is transcribed and translated into gene products x, 
with a delay r. In this example, gene products form 
dimers that act as transcriptional repressors, in¬ 
hibiting transcription of the gene (blunted arrow), 
and decay at a rate c. (B) and (C) Numerical so¬ 
lutions of Eq. © describing the oscillator in A. 
(B) Phase space: monomer concentration x(t) vs. 
the delayed concentration x{t — r) settled to a limit 
cycle. (C) Dimer concentration oscillates as a func¬ 
tion of time. Parameters in B, C: bP = 2, r = 1, 
c = 1, aj 0 = 1, AT = 12, M = 6. 


cooperative binding of transcriptional regulators in 
the case of the zebrafish segmentation clock. There¬ 
fore, although cooperative binding is not ruled out, 
this lack of evidence raises the general question of 
what contribution could be expected from the mul¬ 
tiple binding sites reported in the promoter of seg¬ 
mentation clock cyclic genes. Here we use theory to 
study how multiple binding sites affect nonlinearity 
and biochemical oscillations in a generic description 
of a genetic oscillator. 


II. A promoter with multiple bind¬ 
ing sites 

We first evaluate the regulatory function of a pro¬ 
moter that contains N binding sites for a transcrip¬ 
tional repressor, see Fig. [2] We consider transcrip¬ 
tional repressors although the results in this section 
are more general and would apply to other types of 
transcription factors. We focus on a single tran¬ 
scriptional repressor for the sake of simplicity, and 
assume that all binding sites are identical. That 
is, binding and unbinding of repressors to the dif¬ 
ferent sites occur at the same rates for all sites. 
Moreover, we assume that there is no cooperativity 
in repressor binding. This means that binding and 
unbinding rates are not affected by the presence of 
bound factors to any of the other sites. 
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Figure 2: Schematic representation of a promoter 
with multiple binding sites, numbered from 1 to 
N. Transcriptional repressors (orange squares) 
bind and unbind from binding sites (numbered 
platforms) at the promoter of a gene (light blue 
stretch). (A) and (B) illustrate two equivalent con¬ 
figurations of the promoter with identical number 
of bound transcriptional repressors having the same 
inhibiting strength. 


060012-2 




















Papers in Physics, vol. 6, art. 060012 (2014) / I. M. Lengyel et al. 


The state Vi of the promoter at any time can 
be characterized by the number i of bound factors, 
which goes from 0 to N. With a rate k + , a free 
repressor binds to an empty site at the promoter, 
which steps from state Vi to state Vi+ 1 ; with a rate 
k~, a bound repressor falls off the promoter, which 
steps from state V t to state Vi- 1 : 



U+ 

V N . ( 1 ) 

k~ 


where xo = k /k + is the equilibrium constant for 
binding of factors. Using the constraint Eq. © we 
express Po in terms of P 

p °= p { i+ £)~ N ■ < 5 > 

Equation © and Eq. ([5]) describe the equilibrium 
occupation of the promoter in terms of the concen¬ 
tration x of the transcriptional repressor. 

III. Abrupt inhibition 


Denoting by Pi the promoter occupation probabil¬ 
ities which describe the fraction of time that the 
promoter spends at state Vi in the thermodynamic 
limit j46] . the kinetics of binding and unbinding 
of transcriptional repressors to the binding sites is 
given by the set of ordinary differential equations 


P 0 = —k + NxPo + k Pi 

Pi = —k + (N — i)xPi — k~iPi (2) 

Pk^ (N — i + l)xPi—i + k (i-\- l)P ? ;+i 

P N = —k NPn + xPn —i , 
together with the conservation law 

N 

^ = E P - (3) 

i=0 

where x is the repressor concentration and P is pro¬ 
portional to the number of gene copies. 

We assume here that binding and unbinding 
of repressors to the promoter occur much faster 
than other processes like transcription, translation, 
transport and decay of molecules. This means that 
the promoter occupation probabilities quickly reach 
equilibrium with a given concentration of transcrip¬ 
tional repressors mm- This situation is described 
by Pi = 0 for all i in Eq. ©, and the resulting 
algebraic equations can be solved by induction to 
obtain 


Pi = 




(4) 


In the previous section, we evaluated the kinetics of 
noncooperative binding to a promoter with multi¬ 
ple binding sites. How does the presence of bound 
transcriptional repressors affect the transcription 
rate of the gene downstream of the promoter? In 
general, the strength of inhibition will depend on 
the number of bound repressors, and the regulatory 
function f(x) will have the form 

N 

f(x) = E a i P i > ( 6 ) 

i=0 

where ao = b is the basal transcription rate in the 
absence of bound repressors and a,i is the transcrip¬ 
tion rate in the presence of i bound repressors. Here 
we assume, for the sake of simplicity, that transcrip¬ 
tion proceeds at its basal rate b while there are M 
or less sites occupied by the repressors, and drops 
to zero when the number of occupied sites is larger 
than M, see Fig. [3] We shall consider an alternative 
scenario below. 

In this situation, Eq. © becomes 
M 

f(x) = b y £p i . (7) 

i—0 

Using Eq. © and Eq. ©, the resulting regula¬ 
tory function for a promoter with N binding sites 
is 
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Figure 3: Abrupt inhibition. Full inhibition occurs 
when more than M sites are bound by transcrip¬ 
tional repressors. (A) Transcription rate as a func¬ 
tion of the number of bound transcriptional repres¬ 
sors. (B) Normalized regulatory function, Eq. 0, 
as a function of repressor concentration x, with 
TV = 12 and M = 0,...,11 (dark blue to dark 
red). 


The zebrafish segmentation clock may be an in¬ 
teresting system to evaluate these results. In the 
herl/her7 locus of zebrafish, an estimated number 
of about 12 binding sites has been reported |321I45| . 
The regulatory functions resulting from Eq. © 
for N = 12 are displayed in Fig. 03. Although 
it is clear that inhibition of the promoter for a 
given level of repressors shifts to the right as M 
increases, it is less obvious how the value of M af¬ 
fects the steepness —that is the nonlinearity— of 
the negative feedback. We use Hill functions to 
parametrize the regulatory functions Eq. (0 in a 
more transparent form, see Appendix. Hill func¬ 
tions are parametrized by a Hill coefficient h char¬ 
acterizing the steepness of the curve, and an inhibi¬ 
tion threshold K that describes the concentration 
of repressors that halves the production rate. We 
fit Hill functions to /n.m and obtain an effective 
Hill coefficient h and effective inhibition threshold 
K , for each value of N and M, Fig. 0\,B. Increas¬ 
ing the number of binding sites N while keeping 
M fixed can increase the Hill coefficient. For fixed 
N , increasing M changes the Hill coefficient in a 
nonmonotonic way: there is an optimal value of M 
that maximizes the Hill coefficient and therefore 
nonlinearity. The effective inhibition threshold K 
changes in a simpler form, increasing both with N 
and M. In conclusion, multiple binding sites can 
effectively increase the nonlinearity of the feedback 
via the regulatory function. 

As discussed above, nonlinearity is an essen- 
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Figure 4: Multiple binding sites can increase non¬ 
linearity and enhance oscillations. (A) and (B) Ef¬ 
fective Hill parameters for the regulatory function 
Eq. 0, with bP = 1 and Xo = 1. (A) Ef¬ 
fective Hill coefficient h. (B) Effective inhibition 
threshold K. (C) and (D) Oscillations described 
by Eq. 0 with bP = 2, r = 1, c=l and xo = 1. 
(C) Amplitude of oscillations. (D) Period of oscil¬ 
lations. In C and D, the white region represents the 
nonoscillatory regime, in which the system settles 
to a fixed point. Color bar labels indicate values in 
each panel. 


tial ingredient in a theory of biochemical oscilla¬ 
tions [15]. We therefore ask how multiple binding 
sites in the promoter affect a biochemical oscilla¬ 
tor. We use the regulatory function Eq. 0 in a 
generic model for genetic oscillations. We consider 
a gene that encodes a protein that forms dimers, 
and these dimers can bind to multiple binding sites 
at the promoter to inhibit transcription, Fig. [lj\. 
We introduce an explicit delay r to account for 
transcription, translation, splicing, and other pro¬ 
cesses involved in the assembly of the gene product 
and its dimerization. We assume that dimerization 
is a fast reaction, with a separation of timescales 
from other processes. Therefore, at any time the 
dimer concentration can be approximated by that 
of the monomers squared. The dynamics of the 
product concentration x(t) is given by 
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where b is the basal production rate, c is the decay 
rate of products, P relates to the number of gene 
copies, r is the total delay for product assembly, 
and xq is the product concentration that halves the 
basal transcription rate. The regulatory function 
Eq. © is parametrized by N and M. This genetic 
oscillator is a reduction of the models proposed by 
Lewis m, and other authors - It describes 

the protein concentration, but does not include the 
mRNA; the duration of transcription and transla¬ 
tion are both included in the delay r. Furthermore, 
it does not describe effects present in theories that 
include more than one regulator [3U[5ni[5I], but 
here it is enough to illustrate the effects of multiple 
binding sites in a simpler context. 

We integrate Eq. © numerically and evaluate 
the resulting dynamics by calculating the ampli¬ 
tude and period of oscillations. In all numerical 
simulations, we use the function dde23 from MAT- 
LAB [22] • Scanning the values of N and M, we 
determine whether the system oscillates in steady 
state: when the difference in the maxima over the 
last ten cycles falls below 0.01, the simulation is 
stopped and we record the output. 

The amplitude of oscillations grows with the 
number of binding sites N, Fig. ©d. The range 
where the system oscillates grows with the num¬ 
ber of binding sites N. The amplitude is nonmono¬ 
tonic in M: for a fixed number of binding sites 
TV, there is an optimal value for M that maximizes 
the amplitude of oscillations, Fig0p. The period 
of oscillations grows with N and decreases with M. 
These results show how the change in nonlinear¬ 
ity observed in the regulatory functions is reflected 
in the oscillations as the number of binding sites 
change. 

IV. Gradual inhibition 

There is strong evidence for the products of some 
cyclic genes of the zebrafish segmentation clock 
binding their own promoters and acting as tran¬ 
scriptional repressors [23I127I32I133II35II44] . However, 


we do not have detailed knowledge of how these 
transcriptional repressors affect transcription rates 
when bound to the promoters of cyclic genes. There 
is evidence from transcriptional analysis of the Hesl 
gene in mouse indicating that inhibition is grad¬ 
ual [30] • While the wildtype Hesl promoter con¬ 
taining all three N Box elements that are bound by 
Hesl proteins showed a 30-fold inhibition of tran¬ 
scription in the presence of Hesl, mutations in one, 
two, and three of the N Box elements showed im¬ 
paired inhibition with 14-, 7-, and 2-fold inhibition, 
respectively [30] , 

Motivated by these results, we consider here a 
scenario where additional bound repressors gradu¬ 
ally reduce transcription rate until it drops to zero, 
Fig. ©V. For the sake of simplicity, we assume that 
transcription rate drops linearly as a function of 
bound repressors, from a basal rate b, to zero for 
k + 1 occupied sites 


( b {l-i/{k + l)) if i < k + 1 ( , 

1 \ 0 if*>fc + l, 1 j 

see Fig. ©4. Using this gradual inhibition in Eq. © 
together with Eq. a and a, we obtain regulatory 
functions fN,k{x) 


f N ,k(x) = bP 



fc +1 

x£ 


2=0 




( 11 ) 


see Fig. ©3. As in the previous case, it is clear that 
the effective inhibition threshold shifts to the right 
as k increases, but it is not so clear if the steepness 
of the regulatory function changes and, if so, how. 
Performing fits to Hill functions, we find that the 
nonmonotonic behavior of the effective Hill expo¬ 
nent h is observed again as k increases, Fig. ©\,B. 
Oscillations are similarly affected by noncoopera¬ 
tive binding with gradual inhibition, Fig. ©d,D. 
These results show that the prediction of an opti¬ 
mal value for the number of bound repressors that 
fully inhibits the promoter is robust with respect to 
the details of how multiple bound repressors reduce 
the transcription rate. 
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Figure 5: Gradual inhibition. Binding of multiple 
transcriptional repressors gradually inhibits tran¬ 
scription. (A) Transcription rate as a function 
of the number of bound transcriptional repressors. 
Transcription occurs at the basal rate b in the ab¬ 
sence of bound repressors, and decreases linearly 
to zero for more than k bound repressors. (B) 
Normalized regulatory function Eq. HID as a func¬ 
tion of repressor concentration x , with N = 12 and 
k = 0,..., 11 (dark blue to dark red). 


V. Discussion 

We studied the effects of multiple noncooperative 
binding sites in the promoter of a genetic oscillator. 
We evaluated the behavior of a promoter with mul¬ 
tiple binding sites when binding of transcriptional 
repressors is noncooperative, Fig. [2] We considered 
two different hypotheses for how bound transcrip¬ 
tional repressors affect transcription rates, Figs. [3] 
and[5] In both cases, we calculated how the number 
of binding sites and the number of bound repres¬ 
sors required to produce full inhibition affect the 
nonlinearity of regulatory functions. We showed 
that there is an optimal value of the number of re¬ 
pressors needed to fully inhibit transcription that 
maximizes nonlinearity, Figs. [3JA.B and|S]\B. This 
increased nonlinearity is reflected in the behavior 
of a genetic oscillator controlled by such regulatory 
functions, Figs.|4pD and[6flD. 

Cooperative binding is a well known means to 
increase the nonlinearity of a biological dynamical 
system l2Tll39l . Here we show that the nonlinearity 
of a regulatory function can be increased by multi¬ 
ple binding sites, even if binding is noncooperative. 
This idea may have application in other biological 
control systems as well. Using the same formu¬ 
lation as we did here, one can also describe non¬ 
linearity in the transcriptional activation of gene 
expression, thereby creating effective on-switches 



Figure 6: Multiple binding sites can increase non¬ 
linearity and enhance oscillations in the presence of 
gradual inhibition. (A) and (B) Effective Hill pa¬ 
rameters for the regulatory function }'N,k- Eq. HEED, 
with bP = 1 and xq = 1. (A) Effective Hill co¬ 
efficient h. (B) Effective inhibition threshold K. 
(C) and (D) Oscillations obtained using regulatory 
functions Eq. (fill) with bP = 2, r = 1, c = 1 and 
xo = 1. (C) Amplitude of oscillations. (D) Pe¬ 
riod of oscillations. In C and D the white region 
represents the nonoscillatory regime, in which the 
system settles to a fixed point. Color bar labels 
indicate values in each panel. 


in developmental and physiological regulatory net¬ 
works. A similar effect has been reported in a theo¬ 
retical study of enzymes with multiple phosphory¬ 
lation sites [551155] . It was found that nonessential 
phosphorylation sites give rise to an increase in ef¬ 
fective Hill coefficients, enhancing ultrasensitivity 
in signal transduction [55] . 

Previous work on the segmentation clock ad¬ 
dressed the case of three binding sites [3(1], mo¬ 
tivated by experimental observations of the Hes 7 
mouse promoter m ■ The authors assumed that 
a single bound dimer inhibits transcription com¬ 
pletely, corresponding to the particular case M = 1 
of the theory developed here. They considered 
both noncooperative and cooperative binding, and 
showed that cooperativity increases the effective 
Hill coefficient as expected. In a follow-up 1551 , 
the authors addressed the case of Hesl regulation 
based on the report of four binding sites in the Hesl 
mouse promoter [30] . Again assuming that a sin- 
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gle bound dimer inhibits transcription completely, 
they used the data from the transcriptional analysis 
of the Hesl gene to estimate an effective Hill coef¬ 
ficient for the Mouse Hesl oscillator, obtaining an 
upper bound of about 3. More recently, the effect 
of two binding sites as compared to a single binding 
site was discussed together with differential decay 
of the monomers m- 

Our theory predicts how changing the number 
of binding sites N, and the number M of bound 
repressors that produce full inhibition affects a sin¬ 
gle cell oscillator. Although an experiment that 
changes the value of M may currently be challeng¬ 
ing in a cellular system, the number of binding sites 
N is more amenable to experimental manipulation. 
For example, binding sites could be mutated [30], or 
deleted from the promoter, or they could be inter¬ 
fered with using genome editing strategies such as 
TALEN [5S] or CRISPR [55] to alter or delete spe¬ 
cific binding sites. To assess the effects of these per¬ 
turbations in experiments may also pose some chal¬ 
lenges. Dropping the number of binding sites from 
N = 12 to N = 6 introduces a period change of 
about 2.5%, while amplitude halves over the same 
range, Fig. HP, D. Experiments will require at least 
such precision to reliably detect changes. 

Our results suggest a possible evolutionary mech¬ 
anism to increase nonlinearity in gene regulatory 
systems. In this mechanism, point mutations in the 
promoter that increase the number of binding sites 
for transcription factors may increase the steepness 
of regulatory functions. If the resulting steeper reg¬ 
ulation performs some function better, such muta¬ 
tions would have a good chance to be conserved 
by natural selection. In the case of the segmenta¬ 
tion clock, the amplitude of oscillations could in¬ 
crease with the number of binding sites, possibly 
reducing the signal to noise ratio. Furthermore, 
the range for oscillations would be wider, making 
the oscillatory regime less sensitive to slow extrin¬ 
sic fluctuations of parameter values. Remarkably, 
it may happen that after an increase in N, an in¬ 
crease in M also raises nonlinearity, see Fig. HJA 
This means that weaker repressors would result in 
better oscillations. This evolutionary mechanism 
would provide a simple way to gradually increase 
the nonlinearity of a feedback or other regulatory 
function. 

The theory for the zebrafish regulatory function 
could be refined using the experimentally measured 




Figure 7: Hill functions are characterized by two 
parameters, the Hill coefficient h and the inhibition 
threshold K. (A) Hill functions with K = 5 and 
h = 1 (red), h = 3 (green) and h = 7 (blue). (B) 
Hill functions with h = 3 and K = 3 (red), K = 5 
(green) and K = 7 (blue). 


relative affinities of the binding sites at the herl 
and her7 promoters [32] • Apart from the effects re¬ 
ported here, the number of binding sites may have 
additional roles. For example, it could serve as a 
buffer for fluctuations in gene expression 
augmenting the precision of genetic oscillations. 
This will be the topic of future work. 
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Appendix: Effective Hill functions 

Hill functions are often used to describe the non- 
linearities present in gene regulatory networks un- 
Hill functions are sigmoidal step functions defined 
by 


/h{x) 


1 

1 + {x/Kf h ’ 


( 12 ) 


where the steepnes of the step is characterized by 
the exponent 2 h, and the inhibition threshold K is 
the concentration of repressor that halves the pro¬ 
duction rate, here scaled to unity, Fig. [7] Here 
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Figure 8: Effective Hill functions Eq. m can fit 
regulatory functions of the multiple binding site 
regulatory function, Eq. © for iV = 12, with 
bP = 1 and Xq = 1- (A-C) Examples of fits of 
Hill functions (solid lines) to regulatory functions 
(dashed lines) for (A) M = 1 (green), (B) M = 6 
(light blue) and (C) M = 11 (purple). (D) Effective 
Hill coefficient h as a function of M. (E) Effective 
inhibition threshold K as a function of M. Dots in 
D and E correspond to fits from panels A, B and 
C. 


we include an explicit factor 2 in the exponent to 
account for the dimerization of the transcriptional 
repressors. One advantage of Hill functions is that 
they are very simply parametrized. In the main 
text, we fit Hill functions to the more complex reg¬ 
ulatory functions Eqs. © and ED- Some fits of 
Eq. (0 in the case N = 12 are displayed in Fig. [8] 
for illustration. 
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